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Following the Fukushima Daiichi Nuclear Power Plant (FDNPP) accident of 2011, which deposited radionu- 
clides across Tohoku and northern Kanto, P-ray dose evaluation has been performed to estimate radiation ex- 
posure for small creatures like insects as well as human skin. Using the Monte Carlo radiation transport code 
MCNP-4C, we calculated the P-ray dose for ^^^"'Te, ^^^Te, ^^^I, ^^^Te, ^^"^Cs and ^^'^Cs in air as a function 
of altitude and in soil. These calculations of p-dose rate for each radionuclide were conducted for the condi- 
tions following the FDNPP accident, with ^^^Cs deposition assumed to be 1000 kBq/m^. Beta-ray dose rate 
was found to be ~ 10-fold (resp. 5 -fold) higher than the y-ray dose rate in the soil (resp. on the ground surface) 
at -20 days after deposition, and ~4-fold (resp. 1.7-fold) higher after 6 months or more. For convenience, the 
height dependence of the ratio for 0, 10, 30, 90, 180 and 365 days after deposition was obtained by a fitting 
function. The cumulative 70 jim P-ray dose at 30, 60 and 90 days after deposition was estimated to be 35, 45 
and 53 mGy for the ground surface, and 61, 79 and 92 mGy in the soil, respectively. These results can be used 
to estimate the external P-ray exposure for small creatures as well as for human skin. 

Keywords: Fukushima Daiichi Nuclear Power Plant accident; P-ray dose; radiotellurium; radioiodine; 
radiocesium 



INTRODUCTION 

The nuclear accident at the Fukushima Daiichi Nuclear 
Power Plant (FDNPP) occurred as a consequence of the 
massive earthquake and associated tsunami that struck Japan 
on 1 1 March 201 1. The radionuclides that were released as a 
result were deposited over a wide area from the Tohoku 
region to the northern Kanto region [1]. Because P-ray s are 
not directly related to effective dose, dose evaluations have 
primarily focused on y-ray s. Beta-rays are easily shielded by 
the surfaces of tissue of even a few millimeters thickness, 
and are therefore usually discussed using the 70 jim dose 
equivalent in the context of skin dose [2, 3]. 



Hiyama et al. reported that the first-voltine adults of the 
pale grass blue butterfly (Zizeeria maha) collected in the 
Fukushima area in May 2011 showed relatively mild abnor- 
malities [4]. More severe abnormalities were observed in the 
Fl offspring from the first-voltine females, which were then 
inherited by the F2 generation. They concluded that artificial 
radionuclides from the FDNPP accident caused physiologic- 
al and genetic damage to this species. 

The eggs and larvae of the pale grass blue are 1-2 mm and 
oblate, or 3-25 mm and prolate; therefore, P-ray s effectively 
transfer their energy through the whole volume of the organ- 
ism. Since larvae of the pale grass blue both crawl on the 
ground on warm days and also burrow into the soil on cold 
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days, they are exposed to P-rays from radioactive cesium, 
iodine and tellurium isotopes deposited on the ground. There 
is therefore a possibility that their reported genetic and physio- 
logical abnormalities might have been caused by radiation, 
such as y- and P-rays from the radionuclides released from the 
FDNPP. To estimate the radiation exposure to small creatures 
and human skin, both the P-ray dose in the air as a function of 
altitude, and the P-ray dose in the soil are reported here. 



MATERIALS AND METHODS 

Monte Carlo transport calculation 

The energy deposition in the air and in ground soil by p-rays 
from deposited radionuclides was simulated using a numerical 




1000 1500 
Beta energy (keV) 

Fig. 1. Beta-ray energy spectra for ^^^""Te, ^^^Te, ^^^I, ^^^Te, 
i^^Csand^^^Cs. 



method proposed and vaHdated previously [5]. The energy 
spectra of P-rays emitted from the radionuclides of interest 
('2*"Te, '^''Te, '^'l, ^^^Te, '^^i, ^^^Cs and "^Cs) were deter- 
mined from the P-decay phase space [6] (Fig. 1). Strontium-90, 
which has long half-life (28.6 y), is an important radionuclide 
for the P-ray exposure of fission products. However, on the 
FDNPP accident, the deposition density ratio of ^^Sr/^^^Cs is 
reported to be 1/100-1/1000 [7], therefore ^^Sr is not taken into 
account in this study. The transport of the p-rays was simulated 
with the Monte Carlo N-Particle transport code version 4C 
(MCNP-4C) [8]. The simulation technique was otherwise iden- 
tical to that published previously [2, 3]. 



Calculation geometry 

Computations of radiation dose in the air and in/on the 
ground by p-rays emitted from radionuclides in the ground 
soil were based on a cylindrical structure (diameter 1 m) [2, 
3], as shown in Fig. 2. The bottom of the cylinder was repre- 
sented by a soil layer (thickness 50 cm), where radionuclides 
were distributed uniformly on the surface at a depth of 5 mm 
and 10 mm. Repeating air layers were modeled above the 
ground up to 5 m (the first air layer was 5 cm thick, others 
were 10 cm thick). A mirror boundary condition was used at 
the lateral wall of the cylinder, shown in Fig. 2, to avoid par- 
ticle loss and reduce computation time. By using this ap- 
proach, the cylindrical geometry could be approximated by 
infinite planes of soil and air. The radiation that escaped 
from the lateral wall was treated as entering the cylindrical 
geometry again at the same location on the lateral wall with 
the same energy and angle as when it escaped. The energy 
deposition in the air layers at various heights above the 
ground was computed similarly to obtain the air dose. 




Fig. 2. Monte Carlo input geometry. 
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Dose calculation 

The p-ray dose at the time of deposition of each radionucHde 
Df (Gy/Bq) is calculated from deposited energy SiQi) (J) in 
each air or soil cell, and the energy deposition in the cell is 
divided by each mass m{h) (kg) at height h and the P-ray 
emission rate per decay as follows: 

D%h) = ^ • a,- (/ ='29"'Te,'29Te/3il, 

m{h) (1) 

'32Te/32i/34cs and ^^^Cs). 

The p-ray dose rate D^- (()iGy/h)/Bq) at the time of deposition 
of each radionuclide is given by 

b',{h)=D%h)-ru-n, (2) 

where is the unit conversion factor for seconds to hours 
(3600 s/h) and is the unit conversion factor for Gy to 



Table 1: Half-life, decay constant, and radioactivity ratio to 
^^^Cs for the FDNPP accident [1, 9, 10, 1 1] 



Nuclide 


Half-life 


Decay 
constant (s" 




Radioactivity 
ratio to ^^^Cs 

Assumed From data 




33.6 d 


2.388 X 10" 


-7 


1.0 


1.0-1.08 




33.6 d 


2.388 X 10" 


-7 


0.7 


0.68-0.73 


131j 


8.02 d 


1.000x10" 


-6 


9.2 


7.3-9.2 




3.204 d 


2.504 X 10" 


-6 


8.3 


5.7-8.3 


132ja 


3.204 d 


2.504 X 10" 


-6 


8.3 


5.6-8.3 




2.065 y 


1.064x10" 


-8 


1.0 


0.95-1.03 




30.04 y 


6.096 X 10" 


-7 


1 


1 



^The half-lives are taken to be those of the parent radionuclide. 



|iGy (10^). Time changes of dose rate bs{t^ h) for each radio- 
nuclide are considered in terms of the decay constant Xi of 
each radionuclide. To obtain the p-ray dose at conditions of 
1000 kBq/m^^^^Cs deposition, such that deposition density 
is similar to that in the region between litate Village and 
Fukushima City, the deposition ratio of radionuclide / to 
^^^Cs and the decay constant are considered: 

bs{t,h) = ^b'{h) -ft • exp(-A,0. (3) 

i 

The ratio fi and decay constant Xf for each radionuclide, 
(/ = i2^"^Te, ^29Te, "^I, "^^e, 134^, "'Cs) were 
taken from the literature [1, 9, 10, 1 1], and the values used are 
listed in Table 1 . 

Calculation of 70 ^im dose 

To estimate the skin dose on a small target at the millimeter 
scale, the depth dose was simulated with the ICRU sphere 
[diameter 30 cm; density 1 g cm~^; elemental composition: 
H (10.1 wt%), C (11.1 wt%), O (76.2 wt%), N (2.6 wt%)] 
geometry. The same radionuclides were assumed to exist on 
a tissue-equivalent plastic sphere (ICRU sphere, diameter 
30 cm), after which P-rays were transported [2, 3]. The 
energy deposition accumulated at each depth, and the attenu- 
ation coefficient at a given depth was computed by using the 
ratio between P-ray dose at the depth of interest and that on 
the surface. The time-dependent attenuation factor was calcu- 
lated using the same technique as in Eq. 3. 

RESULTS 

Beta-ray air dose 

The calculation results of the P-ray dose per P emission and 
the P-ray dose rate (jiGy/h) for 1 Bq for each nuclide Z)f as a 
function of height in a condition of 5 mm uniform source 
geometry are shown in Fig. 3a and b, respectively. 
Iodine- 132, ^^^"^Te and ^^^Te, the P-rays of which have 
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Fig. 4. Dose rate change as a function of height from the ground in a condition of 1000 kBq/m Cs deposition, (a-f ) indicate 0, 10, 30, 
90, 180 and 365 days after deposition, respectively. 



relatively high maximum energies (Fig. 2), occur in higher 
doses on the surface and also show height dependence. In 
contrast, low-energy P-rays from ^^^Te show a rapid decrease 
with height. 

The P-ray dose rate changes with time from deposition. To 
show the order of dose rate for the FDNPP accident, the 
p-ray dose rate for 1000 kBq/m^ of ^^^Cs deposition with a 



5 mm uniform source was calculated according to Eq. 3. 
This concentration is similar to that in the region between 
litate Village and Fukushima City. The calculation results of 
the dose rate as a function of height from the ground given 
1000 kBqW of ^^^Cs deposition are shown in Fig. 4a-f, 
which correspond to 0, 10, 30, 90, 180 and 365 days after de- 
position, respectively. 
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Iodine- 132 is the primary contributor to P-ray air dose for 
several days after deposition, followed by ^^^1. These two 
radioactive iodine isotopes have short half-lives (3.204 and 
8.02 days, respectively) as listed in Table 1. Therefore, about 
30 days after deposition the contributions of both ^^^I and ^^^I 
to the dose rate lowered, at which time ^^^"^Te and ^^^Te 
instead constituted the majority of the air dose. After several 
months, the dominant radionuclides were ^^"^Cs and ^^^Cs due 
to their long half-lives (2.06 and 30.02 years, respectively). 

For a 5 mm uniform source, the dose rate in air at the 
ground surface at 0, 10, 30, 90, 180 and 365 days after de- 
position was 240, 81, 33, 16, 11 and 9.7 jiGy/h, respectively. 
This shows that much higher P-dose rates in air than y rays 
were estimated by this calculation. 



Beta-ray dose in soil and on the ground surface 

The calculated P-ray dose rates in soil and on the ground 
surface as a function of time are shown in Fig. 5a and b. 
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Fig. 5. Dose rate change according to the time after deposition in 
the case of 1000 kBq/m^ ^^^Cs deposition (a) in soil, and (b) at 2.5 
cm from the ground surface. 



respectively. The dose rate in soil at deposition is almost 
430 jiGy/h, and this fell to ~ 18 jiGy/h after 1 year. The P-ray 
dose rate in soil was 1.7-1.9-fold higher than that on the 
ground surface. 



Ratio of P-ray dose to y-ray dose 

A P-ray to y-ray dose ratio is useful for reconstructing the 
P-ray dose from the y-ray dose survey monitoring data. The 
y-ray dose value is taken from Imanaka et al. [9]. The P-ray 
to y-ray dose ratio has been obtained for conditions in soil, 
on the ground surface, and 50 and 100 cm above the ground, 
as shown in Fig. 6a. The ratio as a function of height is 
shown in Fig. 6b. The estimated P-ray dose rate is ~ 10-fold 
(resp. 5-fold) higher than the y-ray dose rate in soil (resp. on 
the ground surface) at ~ 20 days after deposition, and ~ 4-fold 
(resp. 1.7-fold) higher after 6 months or more. 

The ratio seen in Fig. 6a has a broad peak at ~ 20 days, but 
after 6 months the ratio levels off. For convenience, the 
height (h) dependence of the ratio rQi) up to 200 cm shown 



(a) 



11 
10 
9 
8 
7 

O 

^ 5 

GO. 

4 
3 
2 
1 
0 





















surface 










— — 50cm 










• 100 cm 












1 











(b) 



50 100 150 200 250 

Time after deposition (day) 



300 



350 



10 




2 0.1 

OQ. 



0.01 



0.001 



■ Oday 

A 10 days 

• 30 days 

□ 90 days 

A 180 days 

O 365 days 



^^^^ 



(3* 



50 



100 
Height (cm) 



150 



200 



Fig. 6. Beta-ray to y-ray dose rate ratio (a) as a function of time 
after deposition, and (b) as a function of height. Curves show the 
fitting results with ^ less than 1 x 10""^. 
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Table 2: Fitted parameters 



Time 

(day)^ 




u 


c 


A 
u 




J 


0 


1.908 


0.0106 


1.917 


0.1608 


0.03427 


-0.0014 


10 


1.864 


0.0112 


2.796 


0.115 


0.03639 


-0.015 


30 


1.533 


0.0119 


2.024 


0.1043 


0.0384 


-0.025 


90 


0.605 


0.0122 


1.538 


0.0959 


0.0389 


-0.0085 


180 


0.189 


0.0140 


1.771 


0.0897 


0.0391 


0.00012 


365 


0.127 


0.0192 


1.687 


0.0896 


0.0401 


0.0012 



Fitting was less than 1x10 ^. ^Days after deposition. 
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Fig. 7. Beta-ray to y-ray dose rate ratio of 10 to 5 mm uniform 
source. 

in Fig. 6b was fitted with parameters of a, b, c, J, e and /as 
follows: 

r{h) aQxp{-bh) + ch~'^Qxp{-eh) +/. (4) 

These fitted parameters are listed in Table 2. The fitted 
curves in Fig. 6b show good reproducibility with the 
obtained p/y ratio. 

Beta-ray air dose ratio of 10 mm to 5 mm 
uniform source 

Calculation of the dose with a 10 mm uniform source was 
also performed, and the result was approximately half the 
value of that calculated with a 5 mm uniform source. The 
ratio between the results of the 10 and 5 mm uniform source 
calculations is shown in Fig. 7. The ratio became slightly 
greater than 0.5 as the height increased, demonstrating that 
high-energy P-rays penetrate to higher locations. However, 
as a height of 200 cm is sufficient for the purposes of dis- 
cussing P-ray dose on insect exoskeleton or human skin, the 
p-ray dose from the 10 mm uniform source could be 



i(r« 




estimated by halving the p-ray dose from the 5 mm uniform 
source. 



Calculation of 70 ^im dose 

The results of P-ray dose per P emission as a function of 
depth for i^^^Te, ^^^Te, ^'^I, ^'^Te, '^^Cs and ^'"Cs are 
shown in Fig. 8. The P-ray dose for each radionuclide is 
decreased by a factor of 10"^ within 1 cm depth, illustrating 
that the P-ray dose of these radionuclides does not contribute 
to the effective dose. 

The attenuation factors as a function of time from 0^00 
days after deposition are shown in Fig. 9a and b. Figure 9a 
shows the attenuation factors for five depths: 70 jam, 450 jam, 
1 mm, 5 mm and 1 cm. Figure 9b shows only the 70 jam at- 
tenuation factor on a linear scale. For convenience, 
least-square fitting was performed for the 70 jim attenuation 
factor with the following equation: 

a{t) = ai + ^2 • exp {(log(r) - asf/al} + 

•exp{(log(0 -^6)V^7}- (5) 

The fitted parameters are listed in Table 3. The 70 jim 
P-ray dose can be estimated by multiplying this attenuation 
factor and the tissue/air kerma ratio (1.09) by the above- 
mentioned surface dose. For example, the cumulative 70 jim 
p-ray dose on the ground surface and in the soil were calcu- 
lated and shown in Fig. 10. The cumulative 70 jim p-ray 
dose in soil is about twice that on the ground surface. It 
shows that insects that burrow into the soil are exposed to a 
higher P-ray dose. The resultant cumulative 70 jim p-ray 
doses at 30, 60 and 90 days after deposition are estimated to 
be 35, 45 and 53 mGy on the ground surface and 61, 79 and 
92 mGy in soil, respectively. 
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Fig. 9. Attenuation factor as a function of time, (a) 70 jam, 450 
jam, 1 mm, 5 mm and 1 cm dose, (b) 70 jam dose and fitting 
function for the data calculated using Eq. 4. 



Table 3: Results of least-square fitting of the 70 (am 
attenuation factor with Eq. 4 



Parameter 


Value 


Error 


ai 


0.220 


0.0392 


a2 


0.0609 


0.0119 


as 


0.5847 


0.0622 


a4 


0.6424 


0.1197 


as 


0.3281 


0.0390 


Cl6 


1.246 


0.0586 


ay 


2.359 


0.2461 




0.00028 





180 




100 150 200 250 300 

Time after deposition (days) 



400 



Fig. 10. Cumulative 70 |im p-ray dose on the ground surface and 
in the soil after deposition in a condition of 1000 kBq/m2 137Cs 
deposition. 



DISCUSSION 

In the case of y rays, ^^"^Cs and ^^^Cs are the main 

contributors to the y-dose rate in air, depending on the time 
after deposition [1, 9]. In contrast, ^^^"^Te and ™Tq also con- 
tribute to the P-ray dose rate in the early stages after depos- 
ition has occurred. The half-lives of ^^^"^Te and ^^^Te (33.6 
d) are slightly longer than those of ^^^I and^^^I (3.204 d and 
8.02 d). When we discuss the P-ray exposure from the depos- 
ited radionuclides of the FDNPP accident, ^^^"^Te and 
™Tq cannot be neglected. On the other hand, ^^Sr, which 
has long half-life (28.6 y), is an important radionuclide for 
the P-ray exposure of fission products. In the case of radionu- 
clides uniformly distributed over a 100 cm^ circular area, the 
skin dose conversion factor for ^^Sr and ^^Y (1.478 and 
1.775 (|iGy-h~^-Bq~^-cm^)) is roughly the same as that for 
^^'^Cs (1.537 (|iGy-h"^-Bq"^-cm^)) [12]. For the FDNPP acci- 
dent, the deposition density ratio of ^^Sr/^^^Cs is reported to 



be 1/100-1/1000 [7] so the skin dose due to ^^Sr-^^Y is 
expected to be 0.1-1% of the skin dose due to ^^^Cs. 
Furthermore, the half-life of ^^Sr (28. 6y) is shorter than that 
of ^^^Cs (30.04y). Therefore, the ^^Sr P-ray contribution is 
negligible in the early stages after the FDNPP accident. 

In this study, the ratio of the P- to y-ray dose rate in air was 
obtained using an empirical fitting function, as indicated in 
Table 2. This result was quite useful for estimating the P-ray 
dose in air from the y-ray dose, because a lot of monitoring 
data has been accumulated for the y-ray dose in air [1, 9, 10, 
13]. However, an assumption about the radioactivity ratios 
has been used in this study (see Table 2); this assumption 
can be used for a relatively wide area, except for part of 
Minami-Soma and Iwaki Cities [1,11]. 

The calculation example of a 70 jim dose (Fig. 10) corre- 
sponds to the skin dose for barefoot walking on 1000 kBq/ 
m^ ^^^Cs-contaminated ground. The resultant cumulative 70 
jim P-ray doses at 30-90 days after deposition are estimated 



Estimation of p-ray dose from Fukushima accident 



483 



to be 35-53 mGy on the ground surface and 61-92 mGy in 
the soil, respectively. While this cumulative P-ray dose is not 
a realistic measure of human exposure, as humans do not 
typically come into barefoot contact with the ground 
throughout the day, it is realistic for insects, which often stay 
in or on the ground. 

CONCLUSION 

This paper has shown the method and results of our P-ray 
dose calculation for the dosimetry of the FDNPP accident. In 
summary, we calculated the p-ray dose in air and in soil for a 
deposition density of 1000 kBqW ^^^Cs after the FDNPP ac- 
cident, using an MCNP simulation. The ^^^""Te, ^^^Te, ^^^I, 
^^^Te, ^^"^Cs and ^^^Cs deposition densities were assumed 
to be in a ratio of 1:0.7:9.2:8.3:8.3:1:1. The P-ray dose rate 
was ~ 10-fold (resp. 5 -fold) higher than the y-ray dose rate in 
soil (resp. on the ground surface) at ~ 20 days after deposition, 
and- 4-fold (resp. 1.7-fold) higher after 6 months or more. 
For convenience, the height dependences of the ratios for 0, 
10, 30, 90, 180 and 365 days were obtained by a fitting func- 
tion, as was the 70 jam attenuation factor. These results can be 
used for p-ray dose estimation, not only for small creatures 
such as insects but also for human skin. The cumulative 70 
\xm P-ray doses at 30, 60 and 90 days after deposition were 
estimated to be 35, 45 and 53 mGy on the ground surface, and 
61, 79 and 92 mGy in the ground soil, respectively. 
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